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Abstract. We present a progress report on a project to derive the evolution of the volumetric 
supernova Type la rate from the Supernova Legacy Survey. Our preliminary estimate of the rate 
evolution divides the sample from Neill et al. [1] into two redshift bins: 0.2 < z < 0.4, and 0.4 < 
z < 0.6. We extend this by adding a bin from the sample analyzed in Sullivan et al. [2] in the range 
0.6 < z < 0.75 from the same time period. We compare the derived trend with previously published 
rates and a supernova Type la production model having two components: one component associated 
closely with star formation and an additional component associated with host galaxy mass. Our 
observed trend is consistent with this model, which predicts a rising SN la rate out to at least z = 2. 
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1. INTRODUCTION 

The importance of Type la supernovae (SNe la) as cosmic distance tracers H, 3] moti- 
vates the efforts to understand their progenitors. We need not only explain the physics 
that makes their explosions so useful, but we must also recognize and control the system- 
atic effects that may result from the properties of their progenitors and their evolution. 



One method used to explore various progenitor models is to compare the evolution 
ofstar formation in the universe with the evolution of the observed SN la rate [H, @, 0, 
ToL 11, 12|) . The challenge of this method is that it requires minimizing systematic 



effects that can produce spurious trends as a function of redshift. The redshift regions 
most important for constraining the progenitors, the highest and lowest redshifts, are also 
the most prone to systematic effects. At low redshifts cosmic variance and the difficulty 
of sampling large volumes and thus produce good statistics are severe challenges. At 
high redshifts, the faintness of the objects produce low signal-to-noise detections and 
less certain typing as well as higher host contamination since projected host offsets are 
smaller. 

The Supernova Legacy Survey (SNLS) offers the opportunity to measure the SN la 
rate trend right in the sweet spot where the volume sampled per square degree is 
large and yet the objects are bright enough for high quality typing. This circumstance 
motivated the recently published study of Neill et al. yj] which produced an anchor point 
for rate evolution studies at z ~ 0.5. Our next step is to extend this effort by binning the 
original sample and adding a higher redshift bin to trace the rate evolution from z = 0.2 
to z = 0.75. The homogeneous sample and high quality typing provided by SNLS can 
thus achieve a reduction in the systematics in this redshift range and provide a foundation 
for rate evolution studies at higher and lower redshifts. 

Even though rigorous constraints of SN la progenitors await accurate low and high 
redshift rate measurements, we can still address a few questions that have arisen from 
other studies of SN la rates. We would like to see if we can reproduce the observed 
trends from previously published studies, in particular those that show a large jump 
in the SN la rate just beyond z = 0.5 fl2|]. Of particular importance to future SN la 
surveys, such as JDEM and LSST, is the question of the rate evolution beyond z = 1, 
which appears to decline dramatically [[Till . This seems to contradict the observation in 
the local universe that the SN la rate is higher in star forming galaxies H, [H, [H, 15] 
which implies that the SN la rate should increase at least out to z ~ 2. This correlation 
of SN la rate with host star formation rate has now been observed in the intermediate 
redshift universe using SNLS [2] which implies it is not a local effect. For this study we 
ask: is the observed trend in the SN la volumetric rate consistent with this correlation of 
rate with host properties? 



2. THE SUPERNOVA LEGACY SURVEY 

The SNLS was instigated with the goal of providing observational constraints on the 
pressure/density ratio of the universe, w. Its first-year results [16], in combination with 



other cosmological probes, provide one of the best observational constraints on w y/Zfl. 
It is designed as a rolling- search survey of four one- square-degree fields evenly spaced 
in RA. Each field is imaged with the Megacam imager on the Canada-France-Hawaii 
Telescope five times a month during its roughly six-month observing season in four 
wavelength bands: g'r'i'z'. A deep u* image is also taken of each field to aid in measur- 
ing host galaxy properties. As each epoch is acquired it is scanned for variable objects 



TABLE 1. Spectroscopically confirmed SN la samples 



Field 0.2<z<0.4 0.4<z<0.6 0.6 < z < 0.75 



Dl 3 13 9 

D2 5 10 12 

D3 6 10 12 

D4 3 8 9 

ALL 17 41 42 



which are recorded in a databaseQ. SN-like variables are assessed for their SN la type 
likelihood II 1811 and all useful SN la candidates are promoted for followup spectroscopy 
with the VLT or the Gemini or Keck Telescopes. The availability of 8 and 10 meter 
class telescopes for spectroscopic followup provides a well defined completeness and a 
large spectroscopically confirmed sample of objects. The database of all variable objects 
provides a means to estimate spectroscopic completeness 01. 



3. RATE EVOLUTION WITH SNLS 

Neill et al. yl used a Monte Carlo method to estimate the SNLS survey efficiency in 



the redshift range 0.2 < z < 0.6. Updated host extinction models H19Q improved the 
estimation of the systematic errors. The study of SN la rates as a function of host 
properties by Sullivan et al. [2] used the same sample and methods, but added an 
additional sample in the range 0.6<z<0.75. 

To estimate the rate evolution using the SNLS, we split the sample from [Q]] and 
[0] into three bins: 0.2 < z < 0.4, 0.4 < z < 0.6, and 0.6 < z < 0.75. Separate Monte 
Carlo runs were performed for each sub-sample including experiments to determine 
the associated systematic errors in each redshift bin. Table 1 lists the spectroscopically 
confirmed sample for each field in each redshift bin. The preliminary results reported 
here will be improved upon by adding more recent SNe la from the survey to increase 
the sample size in each bin. 

We applied the same methods used in yj] to derive the volumetric SN la rate in each 
bin, and associated systematic and statistical errors. We now compare the results to the 
two-component model of SN la production 



4. THE TWO-COMPONENT MODEL 

This model proposes a component associated with the host galaxy mass, to explain the 
non-zero SN la rate observed in early-type galaxies, and a component associated with 
the host galaxy star formation rate, to explain the correlation of SN la rate with host star 
formation rate iS Il5|1 . The volumetric SN la rate evolution, ry(t), is thus expressed as 



see |http : / / legacy . astro . utoronto . ca| and |http : / /makiki . cf ht . hawaii . edu : 872 /sne/| 
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FIGURE 1. Observed SNLS SN la volumetric rate evolution. The filled squares are the SNLS rates 
in each of the three redshift bins. The filled circles are previously published SN la rates derived from 
sam ples primarily confirmed by spectroscopy from the following references (in redshift order): ll3L fsL fToL 
HI &UH 0, [Oh- The open circles are the SN la rates from fTHl with samples having only 50% spectroscopic 
confirmation. The open triangles are the rates from lfl2h . whose samples are confirmed almost entirely with 
photometric methods. Error bars represent the systematic and statistical errors added in quadrature. 



follows: 

rv(t)=AM*(t)+BM*(t), (1) 

with A in terms of SNe la per year per unit host stellar mass and B in terms of SNe la per 
year per unit host star formation. The B component produces SNe la on shorter (~ lGyr) 
time-scales, while the A component is responsible for the longer time-scale required to 
produce SNe la in galaxies dominated by older stellar populations. 

Figure 1 shows a fit (solid line) of the SNLS rate evolution (filled squares) and various 
rates from the literature (solid circles) to the two-component model based on the star 
formation history (SFH) as parameterized by Hopkins and Beacom Q2JJ] . The resulting 
coefficients are annotated on the plot and the resulting evolution of each component 
is indicated by the dashed line (B component) and the dotted line (A component). The 
statistical uncertainty in the combined components is indicated by the hashed region. 

We note that this model appears to be incompatible with two features in the observed 
rates. This fit is statistically inconsistent with the large jump in the rates seen at z = 0.5 
by Barris and Tonry fll2fl. It also predicts an increase in the SN la rate out to z > 2, while 
the rates from Dahlen et al. [11] show a steady decline after a redshift of z ~ 1. The 
highest redshift observation (z = 1 .6) is only marginally inconsistent with the model. 



The mass used in fitting the A and B components was derived by integrating the SFH 
from Hopkins and Beacom 112 ill , and thus includes the mass of dead stars. When fitting 
host properties, the mass derived does not include dead stars. Our fitted values for A and 
B are consistent with published values 01 §6.1.1], but are marginally offset from host 
property values as expected from the differences in host masses. Our final results will 
include a mass correction for dead stars and thus produce A and B values more consistent 
with those derived from host properties. 

5. CONCLUSIONS 

By restricting our measurements of the SN la volumetric rate evolution to samples 
whose members are primarily identified with spectroscopy, we find a trend that appears 
compatible with the observed properties of SN la hosts. The most stringent test of this 
model awaits further measurements of the rate beyond z = 1, but we must be careful. 
If the two-component model is correct, the higher redshift SNe la are more closely 
associated with star formation. This close association could introduce systematics from 
the dust and other features of high-redshift star formation that compound the difficulty 
of deriving accurate high-redshift SN la rates. This caution applies to the use of these 
high-redshift SNe la for cosmology as well. 
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